Deconvolution of dynamic mechanical networks 
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Time-resolved single-molecule biophysical experiments yield data that contain a wealth of dy- 
namic information, in addition to the equilibrium distributions derived from histograms of the time 
series. In typical force spectroscopic setups the molecule is connected via linkers to a read-out de- 
vice, forming a mechanically coupled dynamic network. Deconvolution of equilibrium distributions, 
filtering out the influence of the linkers, is a straightforward and common practice. We have devel- 
oped an analogous dynamic deconvolution theory for the more challenging task of extracting kinetic 
properties of individual components in networks of arbitrary complexity and topology. Our method 
determines the intrinsic linear response functions of a given molecule in the network, describing the 
power spectrum of conformational fluctuations. The practicality of our approach is demonstrated 
for the particular case of a protein linked via DNA handles to two optically trapped beads at con- 
stant stretching force, which we mimic through Brownian dynamics simulations. Each well in the 
protein free energy landscape (corresponding to folded, unfolded, or possibly intermediate states) 
will have its own characteristic equilibrium fluctuations. The associated linear response function is 
rich in physical content, since it depends both on the shape of the well and its diffusivity — a measure 
of the internal friction arising from such processes like the transient breaking and reformation of 
bonds in the protein structure. Starting from the autocorrelation functions of the equilibrium bead 
fluctuations measured in this force clamp setup, we show how an experimentalist can accurately 
extract the state-dependent protein diffusivity using a straightforward two-step procedure. 



Force spectroscopy of single biomolecules relies most 
commonly on atomic force microscope (AFM) [lHZ| or op- 
tical tweezer [8l4l2lll4| techniques. By recording distance 
fluctuations under applied tension, these experiments 
serve as sensitive probes of free energy landscapes 0- 
llH fl4J , and structural transformations associated with 
ligand binding or enzymatic activity [3, @, 0| • All such 
studies share an unavoidable complication: the signal of 
interest is the molecule extension as a function of time, 
but the experimental output signal is an indirect measure 
like the deflection of the AFM cantilever or the positions 
of beads in optical traps. The signal is distorted through 
all elements in the system, which in addition typically 
include polymeric handles such as protein domains or 
double-stranded DNA which connect the biomolecule to 
the cantilever or bead. As shown in the case of an RNA 
hairpin in an optical tweezer [14j . handle fluctuations 
lead to nontrivial distortions in equilibrium properties 
like the energy landscape as well as dynamic quantities 
like folding/unfolding rates. If an accurate estimate of 
the biomolecule properties is the goal, then one needs a 
systematic procedure to subtract the extraneous effects 
and recover the original signal from experimental time 
series data. 

Static deconvolution, which operates on the equilib- 
rium distribution functions of connected objects, is a 
well-known statistical mechanics procedure and has been 
successfully applied to recover the free energy landscape 
of DNA hairpin s [9 , [u| and more recently of a leucine 
zipper domain [14]. In contrast, for dynamic proper- 
ties of the biomolecule, no comprehensive deconvolution 
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method exists. Handles and beads have their own dissi- 
pative characteristics and will tend to suppress the high- 
frequency fluctuations of the biomolecule and as a result 
distort the measured power spectrum. In the context of 
single-molecule pulling experiments, theoretical progress 
has been made in accounting for handle/bead effects on 
the observed unfolding transition rates [la - tlq . How- 
ever, the full intrinsic fluctuation spectrum, as encoded 
in the time-dependent linear response function, has re- 
mained out of reach. The current work presents a sys- 
tematic dynamic deconvolution procedure that fills this 
gap, providing a way to recover the linear response of a 
biomolecule integrated into a mechanical dissipative net- 
work. We work in the constant force ensemble as appro- 
priate for optical force clamp setups with active feedback 
mechanisms [13, [2(| or passive means [12, [21| . While our 
theory is general and applies to mechanical networks of 
arbitrary topology, we illustrate our approach for the spe- 
cific experimental setup of Ref. [14|: a protein attached 
to optically trapped beads through dsDNA handles. The 
only inputs required by our theory are the autocorre- 
lation functions of the bead fluctuations in equilibrium. 
We demonstrate how the results from two different exper- 
imental runs — one with the protein, and one without — 
can be combined to yield the protein linear response func- 
tions. 

We apply this two-step procedure on a force clamp 
setup simulated through Brownian dynamics, and ver- 
ify the accuracy of our dynamic deconvolution method. 
Knowledge of mechanical linear response functions forms 
the basis of understanding viscoelastic material proper- 
ties; the protein case is particularly interesting because 
every folding state, i.e. each well in the free energy land- 
scape, will have its own spectrum of equilibrium fluctua- 



tions, and hence a distinct linear response function. Two 
key properties determine this function: the shape of the 
free energy around the minimum, and the local diffusiv- 
ity. The latter has contributions both from solvent drag 
and the effective roughness of the energy landscape — 
internal friction due to molecular torsional transitions 
and the formation and rupture of bonds between differ- 
ent parts of the peptide chain. The diffusivity profile 
is crucial for getting a comprehensive picture of protein 
folding kinetics and arguably it is just as important as 
the free energy landscape itself for very fast folding pro- 
teins [22U24J . Our dynamic deconvolution theory pro- 
vides a promising route to extract this important protein 
characteristic from future force clamp studies. 



I. RESULTS AND DISCUSSION 



A. Force clamp experiments and static 
deconvolution 



As a representative case, in this paper we consider the 
double trap setup shown in Fig. HJa), which typically in- 
volves two optically-trapped polystyrene beads of radius 
~ O(10 2 nm), two double-stranded DNA handles, each 
~ O(10 2 nm), attached to a protein in the center (l4|. 
For fixed trap positions and sufficiently soft trapping po- 
tentials, the entire system will be in equilibrium at an 
approximately constant tension F. We are interested in 
a force regime (F > 10 pN in the system under consid- 
eration) where the handles are significantly stretched in 
the direction parallel to the applied force (chosen as the 
z axis), and rotational fluctuations of the handle-bead 
contact points are small. Since the experimental setup 
is designed to measure the z separation of the beads as 
a function of time, we focus entirely on the dynamic re- 
sponse of the system along the z direction. However, the 
methods below can be easily generalized to the transverse 
response as well. Though we consider only a passive mea- 
surement system in our analysis, an active feedback loop 
that minimizes force fluctuations can also be incorpo- 
rated, as an additional component with its own charac- 
teristic dynamic response (with the added complication 
that the response of the feedback mechanism would have 
to be independently determined). 

To set the stage for our dynamic deconvolution the- 
ory we first illustrate the static deconvolution for two 
objects X and Y connected in series under constant 
tension F, e.g. a protein and a handle. Let V x {z) 
and V Y (z) be the constant- force probability distribu- 
tions for each of these objects having end-to-end dis- 
tance z. The total system end-to-end distribution is given 
by V XY {z) = jdz'V x (z')V Y (z - z'). In terms of the 
Fourier-transformed distributions, this can be stated sim- 
ply through the ordinary convolution theorem, V (k) — 
V x (k)V Y (k). If 'P XY is derived from histograms of the 
experimental time series, and if V Y can be estimated 
independently (either from an experiment without the 
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FIG. 1. (a) Double optical tweezer force clamp setup for the 
study of equilibrium protein dynamics, with soft traps ap- 
proximating a constant tension F. (b) To define linear re- 
sponse functions, consider an individual component at ten- 
sion F. A small additional oscillatory force /l exp(— iuit) ap- 
plied at the left end leads to endpoint oscillations with am- 
plitudes z L = J 3 eif,L(w)/L and z R — J C ross(^)/L which de- 
fines the self and cross response functions, (c) Two objects 
X and Y connected in series behave as a composite object 
XY whose response functions can be derived through simple 
rules (Eq. |[2|) from the individual X and Y response func- 
tions, (d) Schematic representation of the optical tweezer 
setup consisting of beads (B), double-stranded DNA handles 
(H) and protein (P), with connecting springs. 



protein, or through theory), then we can invert the con- 
volution relation to solve for the protein distribution V x 
and thus extract the folding free energy landscape. A 
similar approach works for multiple objects in series or 
in parallel. 



B. Dynamic response functions 

Before we consider dynamic networks, we define the 
linear response of a single object under constant stretch- 
ing force F along the z direction, as shown in Fig. QJb). 
Imagine applying an additional small oscillatory force 
/l exp(— iuit) along the z axis to the left end. The result 
will be small oscillations z^exp(— iuit) and ZR,exp(— iuit) 
of the two ends around their equilibrium positions. The 
complex amplitudes zl and zr are related to /l through 
linear response: z L = J sc if,L(w)/ L , z R = Jcross(w)/ L , 
defining the self response function J S eif,L(w) of the left 
end and the cross response function J cross (uj). If the 



oscillatory force is applied instead at the right end, 
the response takes the form: zr = J S eif,R(w)/R, zl = 
Jcross(w)/R. Note that since the object is in general 
asymmetric, J se if,L( a; ) an d J S eif,R.(w) are distinct func- 
tions. However, there is only a single cross response (in 
the absence of time-reversal breaking effects such as mag- 
netic fields [25]). For the purposes of dynamic deconvolu- 
tion of a network, these three response functions contain 
the complete dynamical description of a given component 
and are all we need. It is convenient to define the end- 
to-end response function J cc (lo), with zr — zl = J co (ui)f, 
where the oscillatory force / exp(— iu)t) is applied simul- 
taneously to both ends of the object in opposite direc- 
tions. This response turns out to be a linear combination 
of the other functions: J cc = J S elf,L + Jse\f,R - 2J cross . 

As an illustration we take the simplest, non-trivial ex- 
ample: two spheres with different mobilities \x\, and /iR 
connected by a harmonic spring of stiffness k. In wa- 
ter we are typically in the low Reynolds number regime 
and an overdamped dynamical description is appropriate. 
If an oscillating force of amplitude /l is applied to the 
left sphere, its velocity will oscillate with the amplitude 
— iiozh = Hh(fh+k[zB,— zi]), while the velocity amplitude 
of the right sphere is given by — itoz-R = — /iRfc^R — zi]. 
Using the above definitions of the response functions we 
obtain 
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where /i = /^l + A*r- By symmetry J se if,R is the same 
as Jseif.L with subscripts L and R interchanged. The 
end-to-end response J oe has a standard Lorentzian form. 
For more realistic force transducers, such as semiflexible 
polymers, J ee will later be written as a sum of Lorentzians 
reflecting the polymer normal modes. Note that when 
k = 0, and the spheres no longer interact, J se if,L = *A*l/w, 
the standard result for a diffusing sphere, and J cr0 ss = 0, 
as expected, since there is no force transmission from one 
sphere to the other. 

Though all the linear response functions are defined 
in terms of an external oscillatory force, in practice one 
does not need to actually apply such a force to deter- 
mine the functions experimentally. As described in the 
two-step deconvolution procedure below, one can extract 
them from measurements that are far easier to imple- 
ment in the lab, namely by calculating autocorrelation 
functions of equilibrium fluctuations. 



C. Dynamic convolution of networks 

Based on the notion of self and cross response func- 
tions, we now consider the dynamics of composites. We 
explicitly display the convolution formulas for combining 
two objects in series and in parallel; by iteration the re- 
sponse of a network of arbitrary topology and complexity 



can thus be constructed. As shown in Fig. [TJc), assume 
we have two objects X and Y connected by a spring. X is 
described by response functions J* ]f L , J* ]f R , and J„ oss , 
and we have the analogous set for Y. The internal spring 
is added for easy evaluation of the force acting between 
the objects, it is eliminated at the end by sending its 
stiffness to infinity. We would like to know the response 
functions of the composite XY object, J^ x , J^ Y , an d 
Across; where the X and Y labels correspond to left and 
right ends, respectively. The rules (with full derivation 
in the Supplementary Information (SI)) read 
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The rules for connecting two objects in parallel are more 
straightforward and read G* Y = G„ + G^, where a is 
any one of the function categories (self, cross or end-to- 
end), and G denote the inverse response functions. One 
particularly relevant realization for parallel mechanical 
pathways are long-range hydrodynamic coupling effects, 
that experimentally act between beads and polymer han- 
dles in the force clamp setup. We derive the parallel 
rule and show an example hydrodynamic application in 
the SI. For simplicity, however, we will concentrate in 
our analysis on serial connections. To proceed, if we set 
X=H and Y=B, we can obtain the response functions of 
the composite handle-bead (HB) object, J™, if we know 
the response functions of the bead and handle separately. 
Our full system in Fig.QJd) is just the protein sandwiched 
between two HB components (oriented such that the han- 
dle ends of each HB are attached to the protein). The 
total system response functions (denoted by "2HB+P") 
in terms of the individual protein and HB functions re- 
sult by iterating the pair convolution in Eq. (J2J) twice. In 
particular, the end-to-end response J^ BtP is given by: 
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This is a key relation, since we show below that both 
J™ B+P and the three HB response functions J™ H , 
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" Li can be derived from force clamp experi- 
mental data. Hence Eq. ^ allows us to estimate the 
unknown protein response function J^ c . 

We note a striking similarity to the signal processing 
scenario [26j, where the output of a linear time- invariant 
(LTI) network (e.g. an RLC electric circuit) is charac- 
terized through a "transfer function". For such networks, 
combination rules in terms of serial, parallel, and feed- 
back loop motifs exist. The result for J^fix m Eq. © 
can be seen in a similar light: the first term is the self- 
response J^x °f object X which is independent of the 



presence of object Y. The rational function in the second 
term is the "feedback" due to X interacting with Y. As 
expected, if the cross-response J* oss connecting the two 
ends of X is turned off, this feedback disappears. In anal- 
ogy to the transfer function theory for LTI systems, our 
convolution rules form a comprehensive basis to describe 
the response of an arbitrarily complicated network inside 
a force clamp experiment. And like the transfer functions 
which arise out of LTI feedback loops, the convolution of 
interacting components consists of a nonlinear combina- 
tion of the individual response functions. The rational 
functions due to the feedback of mechanical perturba- 
tions across the connected elements are non-trivial, but 
can be exactly accounted for via iteration of the convo- 
lution rules. 



D. Two-step dynamic deconvolution yields protein 
dynamic properties 

To illustrate our theory, we construct a two-step pro- 
cedure to analyze the experimental system in Fig. UJa) , 
with the ultimate goal of determining dynamic protein 
properties in the force clamp. Under a constant force F, 
the protein extension will fluctuate around a mean cor- 
responding to a folded or unfolded state. Though it has 
been demonstrated that under appropriately tuned forces 
the protein can show spontaneous transitions between 
folded and unfolded states, we for the moment neglect 
this more complex scenario. (In the SI we analyze sim- 
ulation results for a protein exhibiting a double- well free 
energy, where the two states can be analyzed indepen- 
dently by pooling data from every visit to a given well; 
the same idea can be readily extended to analyze time 
series data from proteins with one or more intermediate 
states.) We consider the protein dynamics as diffusion of 
the reaction coordinate z P c (the protein end-to-end dis- 
tance) in a free energy landscape C/p(z P e ). The end-to- 
end response function J^ e reflects the shape of Up around 
the local minimum, and the internal protein friction (i.e 
the local mobility), which is the key quantity of inter- 
est. The simplest example is a parabolic well at position 
z p = z Q , namely t/ P (z p ) = Up{z ) + {l/2)kp(z^ c -z ) 2 . If 
we assume the protein mobility /ip is approximately con- 
stant within this state, the end-to-end response is given 
by 
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which has the same Lorentzian form as the harmonic 
two-sphere end-to-end response in Eq. ([T]). Depending 
on the resolution and quality of the experimental data, 
more complex fitting forms may be substituted, includ- 
ing anharmonic corrections and non-constant diffusivity 
profiles (an example of these is given in the two-state 
protein case analyzed in the SI). However for practical 
purposes Eq. Q is a good starting point. 



An experimentalist seeking to determine J^ c (uj) would 
carry out the following two-step procedure: 

First step: Make a preliminary run using a system 
without the protein (just two beads and two handles, 
as illustrated in Fig.^a)). As described in the Materials 
and Methods (MM), time derivatives of autocorrelation 
functions calculated from the bead position time series 
can be Fourier transformed to directly give J s 2 e ™ and 
J 2 ™. The convolution rules in Eq. © relate J®?P and 
J 2 r ™ to the bead/handle response functions, J™ and 
Across! which via another application of Eq. @ are re- 
lated to the response functions of a single bead and a 
single handle. The bead functions J^ if and J^ oss depend 
solely on known experimental parameters (MM) , leaving 
only the handle functions J^ lf and J^ oss as unknowns 
in the convolution equations. Choosing an appropriate 
fitting form, determined by polymer dynamical theory 
(see MM), we can straightforwardly determine J^ lf and 
7 H 

w cross" 

Second step: Make a production run with the protein. 
Eq. © relates the resulting end-to-end response J^ B+F , 
extracted from the experimental data, to the response of 
the protein alone J^ e . Since the first step yielded the com- 
posite handle-bead functions J™ H , </™b> Across which 
appear in Eq. ©, the only unknown is J^ c . We can 
thus solve for the parameters /ip and fcp which appear in 
Eq. g]). 

This two-step procedure can be repeated at different 
applied tensions, revealing how the protein properties 
(i.e. the intramolecular interactions that contribute to 
the diffusivity /ip) depend on force. Even analyzing the 
unfolded state of the protein might yield interesting re- 
sults: certain forces might be strong enough to destroy 
the tertiary structure, but not completely destabilize the 
secondary structure, which could transiently refold and 
affect /ip. 



E. Simulations validate the deconvolution 
technique 

To demonstrate the two-step deconvolution procedure 
in a realistic context, we perform Brownian dynamics 
simulations mimicking a typical force clamp experiment: 
two beads that undergo rotational and translational fluc- 
tuations are trapped in 3D harmonic potentials and con- 
nected to two semiflexible polymers which are linked to- 
gether via a potential function that represents the pro- 
tein folding landscape (see MM for details). We ignore 
hydrodynamic effects which can easily be accounted for 
through parallel coupling pathways, as mentioned above. 

We begin with the first step of the deconvolution pro- 
cedure. A snapshot of the simulation system, two handles 
and two beads without a protein, is shown in Fig. [^a) . 
A representative segment of the z 2 ™(£) time series is 
shown in Fig. [5Jb) . Equilibrium analysis of the time se- 
ries yields the end-to-end distribution 7 ,2HB (z 2 ( J iB ), which 
is useful for extracting static properties of the protein like 
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FIG. 2. The first step in the deconvolution procedure, needed to to determine the handle-bead response J HB : analysis of the 
optical tweezer system without the protein, (a) From top to bottom: a Brownian dynamics simulation snapshot; schematic 
representation of the system; after the first convolution step handles and beads are grouped into composite handle-bead (HB) 
objects; after the second convolution step the full system ("2HB") constitutes a single object, (b) Part of the simulation time 
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series for the total system end-to-end distance z ee . The time series yields the equilbrium probability distribution V (z, 
shown on the right, (c) Top: the MSD functions A 2 ™^) and A™ B (£) (Eq. ([5])) calculated from the simulation; bottom: the 
time-domain response functions J™ B (i) = (/?/2)dA™ B (t)/dt, J 2 e HB (i) = (/3/2)dA 2H:B (£)/dt. Simulation results (symbols) are 
numerical derivatives of the curves in the top panel. The solid lines are a 5-exponential fit to the simulation results, (d) Left 
column: the real and imaginary parts of the J 2HB self and end-to-end response functions. Simulation results (symbols) are just 
the Fourier transforms of the multi-exponential fits in (c). Theoretical fitting results according to Eq. @ and based on HB 
functions J HB are shown as solid lines. Right column: the HB response functions, as determined by the theoretical fitting to 
the full system data. 



the free energy landscape: when the protein is added to 
the system, the total end-to-end distribution is just a 
convolution of the 2HB and protein distributions. (The 
asymmetry of 7? 2HB seen in Fig.[5Jb) arises from the semi- 
flexible nature of the handles.) As described in MM, we 
use the time series to calculate self and end-to-end MSD 
curves A^ f B (£) an d Alf B (t) [Fig.^c)] whose derivatives 
are proportional to the time-domain response functions 
J s 2 o?f B ( i ) and J ec HB W- Tlie multi-exponential fits to these 
functions are illustrated in Fig. [He), and their analytic 
Fourier transforms plotted in the left column of Fig.[2jd). 
We thus have a complete dynamical picture of the 2HB 
system response. However, in order to use Eq. ([3]) to ex- 
tract the protein response, we first have to determine the 
handle-bead response functions J™ . Although a general 



fit of J™ is possible, it is useful to apply the knowledge 
about the bead parameters and symmetry properties of 
the handle response. The handle parameters (the set 
{k B : fi^} in MM Eq. ([7J) are the only unknowns in the 
three HB response functions: J™ H , >/™b> and ^ee B - 
Note that the handle-bead object is clearly asymmetric, 
so the self response will be different at the handle (H) 
and bead (B) end. Convolving two HB objects accord- 
ing to Eq. © and fitting the handle parameters to the 
simulation results for J s 2 e ™(aj) and J™ b (lu) leads to the 
excellent description shown as solid lines on the left in 
Fig. Hid). The handle parameters derived from this fit- 
ting completely describe the HB response functions J™ , 
shown in the right column of Fig. ^d) . The HB response 
curves reflect their individual components: there is a low 
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FIG. 3. Real (top) and imaginary (bottom) parts of the total 



end-to-end response J e , 



(cj) of an optical tweezer system 



with the protein modeled as a single parabolic potential well 
(fep = 0.02 ksT/a 2 , /ip — 0.05/io)- Symbols are simulation 
results, and the solid line is the theoretical prediction, based 
on the convolution of the protein response Jg e (uj) with the HB 
response functions J HB of Fig.[2][d) according to Eq. ([3]). For 
comparison, J^ci^) (dashed line) and J^™(uj) (dot-dashed 
line) are also included. Insets: to show the sensitivity of the 
theoretical fitting, zoomed- in sections of J™ B+F (to) near the 
maxima of the real (top) and imaginary (bottom) compo- 
nents. Both simulation (symbols) and theoretical (blue/red 
curve) results are plotted. The thin pink/cyan curves are the- 
oretical results with /ip different from the true value: from left 
to right, p P = 0.01, 0.03, 0.07, 0.09/i . 



frequency peak/plateau in the imaginary/real part of the 
HB self response, related to the slow relaxation of the 
bead in the trap. The higher frequency contributions are 
due to the handles, and as a result they are more promi- 
nent in the self response of the handle end than of the 
bead end. There is a similarly non-trivial structure in 
the end-to-end response, due to the complex interactions 
between the handle normal modes and the fluctuations 
of the trapped bead (see SI for more details) . 

The double-HB end-to-end distribution in Fig. HJb) 
and the HB response functions in Fig.^d) are all we need 
to know about the optical tweezer system: the equilib- 
rium end-to-end distribution and linear response of any 
object which we now put between the handles can be re- 
constructed. We will illustrate this using a toy model of 
a protein. In our simulations for the second step, we 
use a parabolic potential Up with Up(z) = fcp = 0.02 
ksT/a 2 , and a fixed mobility /ip = 0.05/xo. Here a = 1 
nm and /zo = l/6irr]a, where r\ is the viscosity of wa- 
ter. This leads to the single-Lorentzian response func- 



tion J^ c of Eq. (0J . If this exact theoretical form of J^ e is 
convolved with the HB response functions from the first 
step according to Eq. ©, we get the result in Fig.[3J very 
close agreement with J™ B+F directly derived from the 
simulated time series data. For comparison we also plot 
the separate end-to-end responses of the protein alone 
(P) and the double-HB setup without a protein (2HB). 
As expected J™ B+F differs substantially from both of 
these as correctly predicted by the convolution theory. 
The effect of adding handles and beads to the protein 
is to shift the peak in the imaginary part of the total 
system response to lower frequencies. Additionally, we 
see in J 2BB+F the contributions of the handle and bead 
rotational motions, which are dominant at higher fre- 
quencies. The sensitivity of the theoretical fit is shown 
in the insets of Fig. [31 zooming in on the maxima of 
Re J e 2 e HB+p and Im J 2 HB ^ P , we plot the true theoretical 
prediction (red/blue curves) and results with /ip shifted 
away from the the true value (thin pink/cyan curves) . In 
fact if fcp and /ip are taken as free parameters, numerical 
fitting to the simulation J 2BB+F yields accurate values 
of: /i P = 0.050/io and fc P = 0.0199 k B T/a 2 . Examples of 
successful deconvolution with other values of the intrin- 
sic protein parameters are given in the double-well free 
energy analysis in the SI. 

In practice, any theoretical approach must take into 
consideration instrumental limitations: most signifi- 
cantly, there will be a minimum possible interval t s be- 
tween data collections, related to the time resolution of 
the measuring equipment. The deconvolution theory can 
always be applied in the frequency range up to u> s = l/t s . 
Whatever physical features of any component in the sys- 
tem that fall within this range, can be modeled and ex- 
tracted, without requiring inaccessible knowledge of fluc- 
tuation modes above the frequency cutoff uj s . In the SI, 
we illustrate this directly on the toy protein discussed 
above, coarse-graining the simulation time series to 0.01 
ms intervals, mimicking the equipment resolution used 
in Ref. |14|. The characteristic frequency of the pro- 
tein within the tweezer setup falls within the cutoff, and 
hence our two-step deconvolution procedure can still be 
applied to yield accurate best-fit results for the protein 
parameters. The SI also includes a discussion of other 
experimental artifacts — white noise, drift, and effective 
averaging of the time series on the time scale t s — and 
shows how to adapt the procedure to correct for these 
effects. 



II. CONCLUSION 

Dynamic deconvolution theory allows us to extract the 
response functions of a single component from the over- 
all response of a multicomposite network. The theory is 
most transparently formulated in the frequency domain, 
and provides the means to reverse the filtering influence 
of all elements that are connected to the component of 
interest. From the extracted single-component response 



function, dynamic properties such as the internal mo- 
bility or friction can be directly deduced. At the heart 
of our theory stands the observation that the response 
of any component in the network is completely deter- 
mined by three functions, namely the cross response and 
the two self responses, which are in general different at 
the two ends. The response of any network can be pre- 
dicted by repeated iteration of our convolution formulas 
for serial and parallel connections. Self-similar or more 
complicated network topologies, as occur in visco-elastic 
media, can thus be treated as well. We demonstrate 
the application of our deconvolution theory for a simple 
mechanical network that mimics a double-laser-tweezer 
setup, but the underlying idea is directly analogous to 
the signal processing rules which describe other scalar 
dynamic networks, such as electrical circuits or chem- 
ical reaction pathways in systems biology [27J]. We fi- 
nally point out that dynamic convolution also occurs in 
FRET experiments on proteins where polymeric linkers 
and conformational fluctuations of fluorophores, as well 
as the internal fluorescence dynamics, modify the mea- 
sured dynamic fluctuation spectrum [28|-|30|. The exper- 
imental challenge in the future will thus be to generate 
time-series data for single biomolecules with a sufficient 
frequency range in order to perform an accurate decon- 
volution. For this a careful matching of the relevant time 
and spatial scales of the biomolecule under study and the 
corresponding scales of the measuring device (handles as 
well as beads) is crucial, for which our theory provides 
the necessary guidance. 



III. MATERIALS AND METHODS 

A. Determining the total system response from the 
experimental time series 

A key step in the experimental analysis is to obtain the 
system response functions J se if(w) and ./cross (w) from the 
raw data (which can either be the double handle-bead 
system with or without a protein). This data consists 
of two time series Zb,l(£) and Zb,r(£) for the left/right 
bead positions from which we calculate the mean square 
displacement (MSD) functions 

1 



A sclf (t) = ^((z B ,LW-^B,L(0)) 2 ) 

+ ^b,r.W-^b,r(0)) 2 ), 
A cc {t) = ((z cc (t)-z cc (0)) 2 ), 



(5) 



From the fluctuation-dissipation theorem [25|, the time- 
domain response functions J sc u(t) and J ee {t) are related 
to the derivatives of the MSD functions: J se if(i) = 
{/3/2)dA scli (t)/dt, J cc {t) = {/3/2)dA cc {t)/dt, where 
(3 = 1/ksT. To get the Fourier-space response, 
the time-domain functions can be numerically fit 
to a multi-exponential form, for example J se if(t) — 
J^j Ci exp(— Ait). In our simulation examples typically 
4-5 exponentials are needed for a reasonable fit. Once 
the parameters Cj and Aj are determined, the expression 
can be exactly Fourier-transformed to give the frequency- 
domain response function, J se if(cj) = y\ Ci/ (Aj — iu). 
An analogous procedure is used to obtain J ee (w). The 
cross response follows as J cm ss = -/self - Jee/2. The power 
spectrum associated with a particular type of fluctuation, 
for example the end-to-end spectrum R cc (uj) (defined as 
the the Fourier transform of the autocorrelation), is just 
proportional to the imaginary part of the corresponding 
response function: R cc (u>) = (2fcsT/w)Im J ee (aj). 



B. Bead response functions 

The response functions of the beads in the optical 
traps are the easiest to characterize, since they depend 
on quantities which are all known by the experimen- 
talist: the trap stiffness fctrap, bead radius R, mobility 
/ib = 1/6ttt)R, and rotational mobility /i rot = l/8irr)R 3 . 
Here r\ is the viscosity of water. For each bead the three 
response functions can be defined as described above, 
with the two "endpoints" being the handle-attachment 
point on the bead surface (zg) and the bead center (zb)- 
The latter point is significant because this position is 
what is directly measured by the experiment. For the 
case of large F, where the rotational diffusion of the bead 
is confined to small angles away from the z axis, the re- 
sponse functions are: 



J. 



s, 1l',b( w ) — 'h 



B 

cross 



M = 



Mb 



MBfctrf 



IUJ 



^sclf\s( w ) 



/'B 



flBk 



trap 



2k B Tfi mt R{F)- 
2fj, TOt RF - iu; 



(6) 



The second term in J^ lf s (ui) describes the contribution 
of the bead rotational motion, which has a characteris- 
tic relaxation frequency 2fi rot RF. This term is derived 
in the SI, though it can also be found from an earlier 
theory of rotational Brownian diffusion in uniaxial liquid 
crystals @. 



where z cc (t) = 2b,r(£) — z B,h(t), and we have averaged 
the self MSD of the two endpoints because they are iden- 
tical by symmetry. Calculating the MSD functions is 
equivalent to finding the autocorrelation of the time se- 
ries: for example, if R cc (t) = (z cc (t) z cc (0)) is the end-to- 
end autocorrelation, the MSD A oc (i) is simply given by 

Aee(t) = 2(ii ee (0) - iie e (i)). 



C. Handle response functions 

The double-stranded DNA handles are semiflexible 
polymers whose fluctuation behavior in equilibrium can 
be decomposed into normal modes. We do not need the 
precise details of this decomposition, beyond the fact that 
by symmetry these modes can be grouped into even and 



odd functions of the polymer contour length, and that 
they are related to the linear response of the polymer 
through the fluctuation-dissipation theorem. From these 
assumptions, the handle response functions have the fol- 
lowing generic form (a fuller description can be found in 
the SI): 



■J£f(«) 

'-'cross ( w j 



i/$ 



JVmodc- 






Wmodo / -|\m,.H 



(7) 



for some set of 2iV mo He + 1 parameters {/x„, fc„}. Note 
that since the handles are symmetric objects, the self 
response of each endpoint is the same function J B lf (w). 
The mobilities (j,^ and elastic coefficients fc B encode the 
normal mode characteristics, with the mode relaxation 
times r n = l/(/i B k B ) ordered from largest (n = 1) to 
smallest (n = N mo d e ). The parameter /i^ is the center- 
of-mass mobility of the handle along the force direction. 
Simple scaling expressions for the zeroth and first mode 
parameters in terms of physical polymer parameters as 
well as the connection between the expressions in Eq. (JTJ 
and Eq. (|7|) are given in the SI. (These same expres- 
sions, with a smaller l p , could describe a completely un- 
folded, non-interacting, polypeptide chain at high force.) 
In practice, the high-frequency cutoff N moc { c can be kept 
quite small (i.e. N mo d e = 4) to describe the system over 
the frequency range of interest. 



D. Numerical inversion of convolution equations 

Care must be taken in manipulating Fourier-space re- 
lationships like Eq. ([3]) . Directly inverting such equations 
generally leads to numerical instabilities due to noise and 
singularities. In our case, we can avoid direct inver- 
sion because the forms of the component functions are 
known beforehand (i.e. Eq. ([6j for the beads, Eq. (O 
for the handles, Eq. Q for the protein). Thus when 
we model the response of the double handle-bead system 
with a protein, J™ B+P (u), we end up through Eqs. (|3J) 
and ([3]) with some theoretical function J 2BB+P (ui, K) 
where K is the set of unknown parameters related to 
the components. Since J 2BB ^ P (lu) is known as a func- 
tion of uj from the experimental time series, we find 



K by minimizing the goodness-of-fit function M(K) = 

E. e n [log | J e 2 e HB ^ P ^)| - log I J e 2 e HB - p ( W , K)|] 2 , where Q 
is a logarithmically spaced set of frequencies, up to the 
cutoff frequency lo s = l/t s determined by the time reso- 
lution t s of the measuring equipment. This is equivalent 
to simultaneously fitting the real and imaginary parts of 
our system response on a log-log scale. 

E. Simulations 

In our Brownian dynamics simulations each handle is 
a semiflexible bead-spring chain of 25 beads of radius 
a = 1 nm, every bead having mobility /iq = l/Gitrja. 
The handle persistence length is l p = 20a. The har- 
monic springs used to connect all components together 
(including the beads making up the handles) have stiff- 
ness 7 = 300 ksT/a 2 . The beads have radius R = 50a, 
and the traps have strength fc tra p = 0.00243 hsT/a 2 , 
which corresponds to 0.01 pN/nm. The traps are posi- 
tioned such that the average force in equilibrium F rj 
3 k B T/a = 12.35 pN. 

To capture the essential features of protein dynamics, we 
construct a simple toy model. The protein is character- 
ized by two vectors: a center-of-mass position r P m and 
an end-to-end separation r P c . Both r P m and the trans- 
verse components of r P obey standard Langevin dynam- 
ics with a mobility /i cm = fj,± = 0.12/io- The internal 
dynamics of protein fluctuations is modeled through the 
longitudinal end-to-end component z P c , subject to a po- 
tential [7p(z p e ) and a mobility /ip. The transverse com- 
ponents (x P e , y P e ) feel a harmonic potential with spring 
constant k± = 1.5 ksT/a 2 . 

The simulation dynamics are implemented through a dis- 
cretized free-draining Langevin equation with time step 
At = 3 x 10~ 4 r, where the time unit r = a 2 /kBTfi n . 
Data is collected every 1000 time steps. Typical simula- 
tion times are ~ O(10 10 ) steps, with averages collected 
from ss 20 independent runs for each system considered. 
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Supplementary Material for "Deconvolution of dynamic mechanical networks" 

I. DERIVATION OF DYNAMIC CONVOLUTION EQUATIONS 

Consider a system under tension F consisting of two objects X and Y connected by a spring of stiffness 7, as shown 
in Fig. 1(c) in the main text. The left and right equilibrium endpoint positions are Zxli ^xr respectively for X, and 
Zyl, Zyr for Y. Imagine applying a small additional oscillatory force /Lexp(— iut) at Zxl along the z direction. In 
the long time limit, every endpoint coordinate in the system would exhibit oscillations of the form z a i exp(— iut) with 
some amplitude z a i, a = X, Y, i = L, R. Let us also denote the instantaneous force exerted by the connecting spring 
as /exp(— iut). In terms of / and the four amplitudes z Q i, we can write down five equations based on the definitions 
of the self/cross response functions for X and Y: 

/ = 7( 2 YL ~ Zxr), 

ZXL = ^selfx/L + ^cross/) 

ZXR = J CIOSS h + ^ S clf,R./> (S.l) 

Z YL = — "/selfi/i 

ZYR = - ./ cr o SS /) 

where the u dependence of the response functions is implicit. We would like to relate the X and Y response functions 
to those of the full system. Since the external force perturbation is applied at the X end of the system, by definition 
zxl = Sx/li zyr = ./cross /l- Solving for zxl and zyr from Eq. (jS.lj) . we find the following expressions for the full 
system response functions in the limit 7 — > 00: 

(JX ) 2 jX jY 

tXY _ tX V,' 7 cross; tXY _ '-'cross ^cross /q 9 \ 

-'self.X — J self,L tX , iY ' '"'cross — T x , rY " V°- Z J 

J sclf,R "+" J self,L J self,R "+" J selfX 

If the force perturbation is applied at YR instead of XL, an analogous derivation yields the remaining self response 
function J^ Y : 

tXY _ tY V^crossj /q q\ 

"^self.Y — J self,R Tx T 7Y ' \ D -°) 

J self,R "+" J self,L 

By iterating Eqs. (|S.2p - (|S.3[) . reducing each pair of components into a single composite object, we can readily obtain 
the convolution equations for an arbitrary number of components in series. 

II. DYNAMIC CONVOLUTION WITH PARALLEL PATHWAYS AND LONG-RANGE 

HYDRODYNAMIC INTERACTIONS 

Though the systems described in the main text consist of multiple objects connected in series, one can extend the 
theory to cases where parallel connections are also present. The most significant experimental realization of such 
connections are objects coupled through long-range hydrodynamic interactions. To make our convolution theory 
comprehensive, in this section we will derive the rules for parallel pathways, and then show in particular how they 
can be used to treat hydrodynamics. 

A. Dynamic convolution for two objects in parallel 

Following the notation of SI Sec. U consider two objects X and Y in parallel (in other words, sharing the same 
left and right end-points). For object a = X, Y, let / a j_, exp(— iut) and / a R,exp(— iut) be the forces that need to be 
applied on that object at the left and right ends in order that the end-points exhibit oscillations z QiJ exp(— iut) and 
z a Rexp(— iut). The relationship between the end-point force and oscillation amplitudes is given by: 

2 qL = ^self,L/ QL ' "/cross/aR) ,„ , N 

(S.4) 

ZaR — J crossJah ~r ^self.RjaR- 
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Since the objects are connected in parallel, we know that zxl = ^yl = -^l and zxr — zyr = Zr- Hence we can write 
Eq. (|S.4[) as a matrix equation, 



where: 



tf a f 

*J L OL1 



\ / p \ , v / Jol ja 

2 L el Joih \ Va _ / ^sclf.L ^cross 



(S.5) 



*=[?), fa - [ J r , J <" = /r^ r? ° SS • (S.6) 

Looking at the full system of two objects together, the total forces on the left and right ends required for a particular 
set of oscillation amplitudes zl and zr are additive: f = fx + fy, where f = (/l, /r). If we define the inverse response 
matrix, tf Q = (tf a ) -1 , then 

f Q = t? Q z (S.7) 

and we obtain: 

f=(^ x + ^ Y )z^ XY z. (S.8) 

This is the main result for convolving parallel pathways: the inverse response matrices of each pathway additively 
combine to yield the total inverse response matrix. Let us label the components of the inverse response matrix G a as 
follows: 

/ -^ ... J3 °"- R ,. n, -^ „. J "°",,_ ™\ /rrc net \ 

(S.9) 

Then since J XY — ( G XY ) -1 , we can solve for the XY response functions in terms of the components of G XY = 

^ x + ^ Y : 

nXY r<XY 

"self.R. .YV u n .„„, 




tXY _ ^self.R ,-XY _ -'.■!■■ .-,.-, 

"self,L ^XY /iXY //nYV \9' "cross 



/oXY nXY _ CriXY \2 ' "cross ^XY /-tXY _ fr»XY \2 

^self.L^sclf.R ^crossJ ^sclf.L^sclf.R Across/ 



r XY ( S - 10 ) 

tXY ^sclf.L 

J self,R ^.XY /-"XY _ (nXY \2 ' 
^self.L^self.R Across J 

Eqs. (|S.9[) - (|S.10[) constitute the complete rules for dynamic convolution of two objects in parallel. 

B. Incorporating hydrodynamics as a parallel connection 

As a simple, realistic application of the above parallel convolution rules, consider the following setup, which consti- 
tutes two parallel pathways: a) two beads of self- mobility /^b trapped in optical tweezers of strength fet rap , interacting 
through long-range hydrodynamics; b) a DNA handle which connects the two beads. (The handle can actually be any 
composite object in series, for example two DNA strands on either side of a protein, like in the main text.) Though 
in principle all the objects in the system are coupled hydrodynamically to each other, we will consider only a single 
pairwise interaction: between the two beads. Since the beads are the objects with the largest drag, this is by far the 
most important interaction for the practical purposes of analyzing experimental data. Though the rotational degrees 
of freedom for the beads can be incorporated in the theory, we will ignore them in this example, since they give only 
a small correction at large forces. Thus the center-of-mass of bead i, and the handle end-point attached to that bead, 
will have the same oscillation amplitude: z%, i = L, R. 

To get the total system behavior, we start by looking at each pathway independently, and calculating the associated 
inverse response matrix. If we assume the handle is symmetric with end-point response functions J^ lf and </„ oss , the 

matrix G H for the handle pathway follows from Eq. (|S.9[) : 



■/,, 



^H = VSSttF (./«,)M^ s ( &J* ^cross . (S . n) 

I "ess -'self / \ Lr_-„„_ (_T = „if / 



(-O^-^SosJ 2 (•0) 2 -(./c H I o B J 2 






12 

For the bead pathway, we can derive G B starting from the equations of motion of the two trapped beads in the 
presence of external force amplitudes /bl and /br acting on the left and right beads: 

-lUZ^ = MB.selfJBL + MB.cioss/bRi ,„ , 

/•tot j-tot (b.IZj 

— JWZr — MB.sclt/BR + MB.cross/BL • 

Here / B °* is the iotaZ force amplitude on the ith bead, including the contribution of the trap force and the external 
force: / B °* = — k tYap Zi + f^. The terms HB,se\t and hb, cross are respectively the diagonal and off-diagonal components 
of the two-bead hydrodynamic mobility tensor. /iB,seif is the bead self-mobility, and /ib, cross is the cross mobility, 
describing the long-range hydrodynamic coupling between the beads. We assume both of these mobilities are roughly 
constant. (The justification for this approximation is given at the end of this section, along with full expressions 
for //B.seif and //b, cross in a typical experimental setup.) We can solve Eq. (JS.12J) for /bj, expressing it in the form 
fB = G B z, where: 



,, /<k ^MB.sclf i<^B,cross \ / „ B •".B 

ttB _ Ktl ' a P (^B, s olf) 2 -(MB.„o s5 ) 2 (MB. s olf) 2 -(MB,„o ss ) 2 1 = I G sclf G d 

^ I ^MB, cross L ^PB.aclf I _ I ffi f<B 

\ (MB, S clf) 2 -(pB,cros S ) 2 Ktl ' a P (MB, S clf) 2 -(MB,„o SS ) 2 / V^Cl'OSS ^sdf 



cross 



(S.13) 



The full system inverse response matrix is t^ HB = ^ H + t^ B . Starting with Eq. (fSTTj) for t^ H and Eq. (|ST3)) for 
C? B , we can derive the components of the total response J HB from Eq. (|S.10[) . 

To complete the description of this system, we need expressions for /ZB,seif and /zb, cross- In equilibrium, the beads 
have some average center-to-center separation AZ, and for added realism, we will consider the bead centers to be a 
height H above a surface (i.e. the microscope slide). In a typical setup AZ and H may be on the order of the bead 
radius R, resulting in non-negligible modifications to each bead's self-mobility from both the presence of the wall 
(assumed to be a no-slip boundary), and the presence of the other bead. Though the mobilities will instantaneously 
depend on the exact positions of the beads relative to themselves and the wall, we will use the mobility values at 
the equilibrium positions, since the amplitude of the bead fluctuations in the traps is small compared to AZ and H. 
Hence the mobilities will be functions of AZ, H , and R. The bead self mobility can be estimated as: 

9R R 3 \ ( 15i? 4 
\m + 8H3J V ~ ~4AZ 4 



MB jS elf = MB(l-:^ + ^) (l-TT^l)> ( S - 14 ) 



where /iB = l/(6irr)R) is the self- mobility of a bead alone in an unbounded fluid. The first parenthesis is the correction 
due to the wall, and the second parenthesis is the correction due to the other bead [1|. Finally, for /zb, cross we use 
the corresponding component of the Blake tensor (at the Rotne-Prager level) , which describes the cross mobility 
of two beads moving parallel to a no-slip boundary: 



_ ^bR 

MB, cross — "Z 



3 



, /16ff 4 - 22H 2 AZ 2 + AZ 
2R Z ' 



V (4H 2 + AZ 2 ) 7 / 2 AZ\ 

(S.15) 



12H 4 + AH 2 AZ 2 + AZ 4 



(A 

\AZ (4H 2 + AZ 2 ) 5 / 2 



III. DETAILS OF THE BEAD RESPONSE FUNCTIONS 

In the absence of bead rotation (where the bead only has center-of-mass translational degrees of freedom), the 
self response functions of both the bead center (zc) and and handle-attachment point on the bead surface (zg) 
are identical: they describe a sphere with translational mobility fiB in an optical trap of strength fc trap , namely 
J B lf B (ui) — J B lf g(w) = hb I Xr^B^trap — iw). (Here, as in the main text, we ignore self-mobility corrections due to the 
microscope slide surface or the presence of the other bead; if desired, these can be simply incorporated as described 
in the previous section.) Similarly, since the bead is a rigid body and any force perturbation at the bead center is 
directly communicated to the bead surface, J\\ oss {oj) is equal to J B lf B ( w ) 



The situation gets more complicated when the rotational degrees of freedom are included, but this involves only the 
surface response function J B | f s , which now has to include an extra term to account for the rotation of zg, along with 
the translational motion of the bead itself. To derive the rotational response (i.e. the second term of J B lf s in Eq. [6] 
of the main text) we take advantage of the fluctuation-dissipation theorem. If R is the vector connecting the bead 
center to the handle-attachment point, such that R z = zs — Zb, and u(i) = ~R(t)/R the corresponding unit vector, 
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then we would like to calculate the time correlation function C R (t) = (R z (t)R z (0)) — (R z ) 2 = R 2 ((u z (t)u z (0)) — (u z ) 2 ) 
in the case where z$ is subject to a constant force F along the z axis. Here (R z ) refers to the equilibrium average 
of R z . For the purposes of analyzing the optical tweezer experiment, we consider only the large force case, where 
F 3> k B T / R, and thus u z (t) is in the vicinity of 1. The fluctuation-dissipation theorem states that the time-domain 
rotational response function is equal to —/3dC R (t)/dt, where (3 = l/fc^T, which we can Fourier transform to get the 
frequency-domain response. 

The key to calculating C R (t) is to note that the equilibrium fluctuations of u(t) correspond to diffusion on the surface 
of a unit sphere subject to the external potential V(u) = —FRu z . Let us define the Green's function G(u, u';i) 
as the conditional probability of ending up at u at time t, given an initial position u' at t = 0, or equivalently 
G(u, u'; 0) = <5(u — u'). This Green's function satisfies the Smoluchowski equation for rotational diffusion Q: 



d_ 
at 



G(u, u'; t) = ix Iot n ■ (KG(u, u'; t) + k B TG(u, u'; t)KV(u)) , 



where Keux d/du. The correlation (u z (t)u z (0)) can be expressed in terms of G as: 



(u,(t)u,(0)) = J d 2 uj d 2 W u z u' z G(u, u'; i)* cq (u'), 



(S.16) 



(S.17) 



where ^I/ eq (u) = Aexp(~/3V(u)) is the equilibrium probability distribution of u, with normalization constant A. 
Taking the time derivative of both sides of Eq. (JS.17J) . and substituting the right-hand side of Eq. (|S.16|) for dG/dt, 
one can derive the following relation describing the time evolution of C R {t) [4J|: 



s c »<'> 



-2k B Tfi mt C R (t) - -^ m tFD R (t). 



(S.18) 



Here D R (t) is the higher order correlation function D R (t) — R 3 ((P2(u z (t))u z (0)) — (P2(u z ))(u z )). where P2(x) = 
(3x 2 — l)/2 is the second-order Legendre polynomial. In fact, dD R (t)/dt can itself be expressed in terms of even 
higher order correlation functions, and this procedure can be iterated to form an infinite set of linear differential 
equations. Rather than formally solving this set, which is possible but tedious, we employ the effective relaxation 
time approximation [4], where C R (t) is assumed to be dominated by a single exponential decay with relaxation time 
r, namely Cn(t) f=a (7fj(0)exp(— t/r). The coefficient Cr(0) = R 2 {(u 2 ) — (u z ) 2 ), and can be evaluated using the 
equilibrium distribution \l/ eq , yielding Cr(0) = (fc^T) 2 jF 2 for large F. Using Eq. (|S.18[) , the relaxation time can be 
written as: 



T = — 



%(0) = 



Cb(0) dt 



c R (p) 



2k B T l i mt C R (0) + -ti mt FD R (0) 



(S.19) 



Like Cij(0), the D R (Q) appearing on the right of Eq. (|S.19|) is an equilibrium average determined through ^ cq : 
Dr(0) = 3(kBT) 2 R/F 2 in the large F limit. Since in this limit the first term in the brackets in Eq. (|S. 19[) is negligible 
compared to the second, we ultimately find t _1 = 2fi rot RF. (The same approximate expression for the relaxation 
time t can also be derived from a theory describing rotational Brownian diffusion in uniaxial liquid crystals (5||.) The 
rotational response is the Fourier transform of —f3dC R /dt = (3C 'fl(0)r _1 exp(— t/r): 



PC R (0)r- 



2k B T^ ot R/F 
2(j, mt RF — iu) 



(S.20) 



IV. DETAILS OF THE HANDLE RESPONSE FUNCTIONS 

We can use standard ideas from the theory of polymer dynamics [3] to derive a simple fitting form for the handle 
response functions. Consider a handle in isolation at equilibrium under constant tension F, described by a continuous 
space curve with z component z(s, t) at time t and contour coordinate s, where s runs from to L. Without deriving 
a detailed dynamical theory, one can still make a few generic assumptions about the behavior of z(s, t). The first is 
that it can be decomposed into a sum over normal modes ^ n (s) in the following way: 



f(»,t)=«d(«,*) + 53Pn(*)*n(*). 



(S.21) 



n=l 



Here the first term zo(s,i) is some reference contour, and the second term represents deviations from that reference 
contour, where P n (t) is the coefficient of the nth normal mode. For convenience we set the reference contour zq(s, t) = 
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(z(s,t)), the thermodynamic average over all polymer configurations in equilibrium. The t dependence of zo(s,i) 
incorporates the center-of-mass motion of the polymer, so that {(zo{s,i) — z o (s,0)) 2 ) = 2D cm t, where D cm is the 
center-of-mass diffusion constant. With this choice of reference contour, (P n (t)) = 0. 

The second assumption is that the equilibrium time correlations of the normal mode coefficients have the form, 

(P„(t)Pm(0)> = S n<m A n exp(-A„t), (S.22) 

for t > 0, where A n is some constant, and A„ is the inverse relaxation time of the nth normal mode. This type of 
exponentially decaying normal mode correlation appears in dynamical theories for many types of polymers: flexible 
chains Q, nearly rigid rods [6|, mean-field approximations for semiflexible chains @, [8|. By convention, we order the 
normal modes such that A n increases with n, i.e. n = 1 corresponds to the largest relaxation time. The third and 
final assumption is that due to the symmetry of the handle (the two end-points are equivalent), the normal modes 
\l/„(s) can be grouped into even and odd functions of s around the center L/2, such that & n (L) = (— 1)"W„(0). 

Putting all these properties together, we can derive expressions for the MSD of a handle end-point, A^ lf (£) = 
((z(L,i)-z(L,0)) 2 ), and the cross MSD, Af IOSS (t) = {(z{L, t) - z(0, 0)) 2 ): 

OO 

Af eU (t) = 2D cm t + 2 £ A n (l - er^l (L), 

n = 1 (S.23) 

A" ro88 (*) = 2D cm t + 2 f>l)"A„(l - e- A "*)* 2 (L). 

n=l 

From the fluctuation-dissipation theorem, the time-domain response functions are related to the MSDs through: 
Ja(t) = (/3/2)dA a (t)/dt, a — self, cross. Taking the Fourier transforms of J^(t) yields the handle fitting forms shown 
in Eq. [7] of the main text: 



n_1 (S.24) 

,-,,H °° I -|Vi,,H V ' 

7 H {, ,\ - !^<L , V^ ^ i ^ En 

crossl ' u ^a^k^-i^ 

n=l r " n 

where u^ = (3D cm , fj^ = /3A n A n ^I' 2 1 (L), k^ — A„//i„, n > 0. In calculating the Fourier transforms, we note that all 
terms in Eq. (|S.23[) are implicitly multiplied by the unit step function 0(£), since the MSD functions are defined only 
for t > 0. 

As a demonstration of the handle fitting forms, consider the simplest case, where a handle consists of two spheres, 
like in the response function example preceding Eq. [1] in the main text. The functions in Eq. [1] can be rewritten in 
terms of a center-of-mass and one normal mode contribution, as expected for a single-spring system: 

■>self,L — 1 ; — , Jself.R — 



/ifc — IU1 U! jjLK — IU! 

Across — ', ; • ^O.ZOJ 

U! /J,K — tU! 

In the symmetric case where /zl = MR = M/2, these response functions have exactly the same form as Eq. IS. 241 with 
(4 = /if = /i/4, fcf = Ak. 

For polymer handles used in experiments there will be contributions from a spectrum of normal modes. However 
the response functions are still dominated by the center-of-mass and lowest-frequency normal mode terms. Hence for 
the purpose of estimating the linear response characteristics of a given setup, we provide simple scaling expressions 
for the parameters u^, uf , and fef in the case of a handle which is a semiflexible polymer of contour length L and 
persistence length l p (i.e. double-stranded DNA). 

By analogy to the two-sphere example above, /1q = /if is just the center-of-mass mobility of the handle along the 
stretching direction, and fcf = 4fc, where k is the effective spring constant of the handle. For large forces, where the 
handle is near maximum extension, we can approximate it as thin rod, for which the mobility /in along the long axis 
is given by /i|| = ln(L / ' d) / '(2irrjL) to leading order. Here 77 is the viscosity of water and d = 2a the diameter of the 
handle. Hence /j^ = /iil~ /tii. To get the effective spring constant, the starting point is the approximate Marko-Siggia 
interpolation formula [9| relating the tension F felt by the handle to its average end-to-end extension z ee along the z 
axis: 



k B T 

F = J h 



4(1 - z oe /LY 



(S.26) 
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Since the force magnitude F is related to the polymer free energy F through F 
spring constant k — d 2 F/dz 2 c — dF/dz cc , we can esth 



dF/dz cc , and the effective handle 



k B T 



1 



1 



L 2L(1 - z ce /Lf 



k fromEq. (|S.26[): 




k B T 
IpL 


V+4 l f - 

\k B T 


L + 4 


4k B T 


fl p F\ 3/2 
\k B Tj ' 




l p L 





3/2 



(S.27) 



where in the second line we have used the fact that F ^$> k B T/l p for the cases we consider, i.e. k B T/l p « 0.1 pN 
for double-stranded DNA, where l p w 50 nm. Since k^ — 4k, we can derive the following expression for the longest 
relaxation time t\ — (/i^fcH)-! of the handle end-to-end fluctuations along the force direction (in other words the 
relaxation time of the first normal mode) : 



T\ = 



irr/lpL 2 



k B T\n(L/d) \k B T 



l P F 



-■ill 



(S.28) 



Up to a constant prefactor, this is the large force limiting case of an earlier scaling expression for the longitudinal 
relaxation time that has been verified through optical tweezer experiments on single DNA molecules [10] . 

With the parameter values a = 1 nm, T — 298 K, and r\ — 0.89 mPa-s, we can estimate a typical relaxation time for 
DNA handles where l p = 50 nm and L ~ O(100 nm): n ~ O(10 ns). For the specific case of the Brownian dynamics 
simulations, where the handles are somewhat shorter at L — 50a, the scaling argument predicts t\ rj 0.7r = 3 ns, 
comparable to the numerically fitted value t± ~ 0.25r — 1 ns, where r = a 2 jk B T [l§ = 4 ns. In Fig. 2(d) in the main 



text, the frequency scale t x s=s 4t 



is where we see the clear divergence between the handle-end and 



bead-end HB self response functions J™^ 1 ^) an d ^™b( w )- This is related to the fact that the handle fluctuations 
become the dominant contribution for the HB object above this frequency scale. (Given the bead radius R = 50a and 
trap strength /ctrap = 0.00243 k B T/a 2 , the characteristic bead frequency scale n B k tra p «5x 10 _5 t _1 = 12500 s _1 is 
much lower.) As a general design principle, experiments will have the best signal-to-noise characteristics when there 
is a good separation of scales between the handle, protein, and bead characteristic frequencies. The example in the 



main text satisfies this idea, since the protein characteristic frequency /ipfcp 
from the ranges of the beads and handles. 



io- 



2.5 x 10 5 s Ms distinct 



DETAILS OF THE INDIVIDUAL HANDLE AND BEAD CONTRIBUTIONS TO THE 
HANDLE-BEAD RESPONSE FUNCTIONS 



As described in the main text, the HB response functions reflect the contributions of their handle and bead com- 
ponents. To elucidate this, we plot in Fig. IS. II the HB functions from Fig. 2, together with the individual H and B 
functions. As expected, the HB self response at the H end is very similar to the individual H response, and at the 
B end it is close to the individual B response. For the HB end-to-end function, only comparison with the handle 
component is possible (since the bead end-to-end response is zero). Again the two functions are similar, but clearly 
adding the bead onto the handle perturbs its end-to-end response, shifting it to lower frequencies and changing its 
shape. The way in which the components contribute to the response of the total object is not merely additive, and 
requires the convolution rules in order to be to accurately predicted. 



VI. DYNAMIC DECONVOLUTION OF A DOUBLE- WELL PROTEIN 



Instead of the single- well protein of the main text, we start with a double- well potential Up(z^ c ) shown in Fig. IS. 21^ ) 
This potential arises from the force: 



+fcizi 



F P (z) 




+ k 2 z 2 



(S.29) 



where the parameters k\ = 2 k B T/a 2 , k m = 1.5 k B T/a 2 , k 2 = 0.15 k B T /a 2 , z\ = 5a, z m = 12.3a, z 2 = —1.3a. The 
energy profile plotted in Fig. IS.2f a) has been additionally tilted by a contribution Fz representing the equilibrium 
stretching at constant force F = 3 k B T / a. To mimic the effects of internal friction on the protein dynamics, we allow 
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FIG. S.l. The panels show the handle-bead (HB) response functions, taken from Fig. 2 in the main text, along with functions 
for the handle (H) and bead (B) separately. The self and end-to-end response functions are in the left and right columns 
respectively, with real parts in the top row and imaginary parts in the bottom row. 

for a coordinate-dependent mobility /ip(zf e ), which is plotted in Fig. IS. 2T b). In the left well /xp = 0.03 /io (folded 
state), while in the right well we have a higher value [ip = 0.12 /i (unfolded state), with a linear transition of width 
2a around the energy barrier. For real proteins, larger internal friction in the folded state may occur due to a greater 
proportion of intact native bonds compared to unfolded configurations. 

Formally, the protein end-to-end response Jq B (uj) may be derived for this double- well potential using a numerical so- 
lution of the corresponding Fokker-Planck equation, based on the method of Bicout and Szabo |lll | . The numerical pro- 
cedure consists of the following steps. We discretize our potential Up(zi) over M values z i} i = 1, . . . , M, in some range 
z-mm to z max . The bin size d = (z max — z m i n )/M. Define local equilibrium probabilities pi = exp(—Up(zi)/kBT)/Z, 
where Z = ^2 t exp(— Up(zi)/k,BT). Transition frequencies w(i\j) of going from bin i to bin j are given by: 



w(i\j) = k B T 



fip(zj) + hp{zj) 
2d 2 



exp 



P(Up{ Zi )~Up(z )) 



(S.30) 



Reflective boundary conditions at the ends of our range are imposed by defining w(0|l) = w(l|0) = w(M\M + 1) = 
w(M + 1\M) = 0. The M X M symmetric rate matrix R has the nonzero elements: 



i? Vi = -w(i + l|i) - w(i - l|i), Ri,i±i = Ri±i,i = \/w(i± l\i)w(i\i± 1) 



(S.31) 



The eigenvalues and eigenvectors of this matrix satisfy an equation of the form: Rv a = — A Q v Q , where the index 
a = 0, . . . , M — 1 and we arrange the eigenvalues in increasing order. With reflective boundary conditions at either 
end, the smallest eigenvalue is always Aq = 0. The correlation function C^ c (t) = (z^ e (t)z^ e (0)) for the protein is then 



given by C c p c (i) - Er=i C« 



-A„t 



where the coefficients C a are: 



C Q — 



M 



fPi&iucti) 



(S.32) 



Here v a i is the ith component of the eigenvector v a . Finally, from the fluctuation-dissipation theorem, the time- 
domain response function J^ c (t) is related to Cf e (t) through J^ c (t) — —(3dC^ c (t)/dt. After a Fourier transform, we 
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FIG. S.2. Analysis of an optical tweezer system with the protein modeled using a double-well potential, (a) The protein 
potential Up as a function of end-to-end distance z co , from Eq. (|S.29|) . The potential as shown is tilted by the external 
stretching force F — 3 knT/a. (b) The protein mobility profile /Up (2^). (c) Formal linear response result for Je e (u)) based on a 
numerical Fokker-Planck solution using the full profiles in (a) and (b). Inset: ImJ c ^(w) (red curve), with analytical estimates 
for the three contributions to the protein response (orange). The sum of these analytical estimates is shown as a dashed black 



, with the portions assigned to either the folded or unfolded wells colored red 

p ) calculated from the 



curve, (d) Fragment of time series for z c 

and blue respectively. The remainder is shaded gray, (e) The full system free energy profile U(z'i 

simulation (dark purple), together with the protein-only result U(z se ) (light purple) obtained by static deconvolution using 

the 2HB distribution in Fig. 2(b) of the main text. This agrees with the original theoretical form, Eq. ([S.29I) . shown as a 

dashed line, (f) For the time series in each well, the corresponding full system potential Ui(z™ B+F ) (blue/red lines) and the 

deconvolved result Ui(z F e ) (cyan/orange lines), i = 1, 2. 



get the response function: 



-£M 



M— 1 Q rt A 
a=l 



A, 



iu> 



(S.33) 



where all the parameters {C Q , A Q } are known numerically. 

This response function is shown in Fig. IS.2f c). Its structure can actually be decomposed into three contributions, as 
shown in the inset: two correspond to fluctuations of the protein about the local minima; the third, at lower frequencies, 
corresponds to transitions between the wells. The orange curves in the inset are simple analytical predictions for these 
contributions: the two intrawell peaks at higher frequencies are just single Lorentzians with the appropriate values of 
k and /Up for each well, and the interwell transition peak is a based on a discrete two-state description [12]. Though 
approximate, the analytical sum (dashed line) captures well the exact numerical result (red line). 

However if we go beyond this and try to naively apply the convolution theory on the Jf c (oj) of Fig. IS.2f c). the results 
for the total system will deviate from the actual behavior as seen in simulations. This is not surprising, since the 
(large amplitude, low frequency) transitions between wells are a fundamentally nonlinear process, and thus violate the 
linear response assumption on which our theory is based. To correctly apply our analysis in this situation, we should 
focus on the linear response of the system within each state, where the protein fluctuates about the local minimum: 
in essence, we break the problem into two single- well systems. 
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FIG. S.3. Analogous to Fig. 3 in the main text, showing the accuracy of the theoretical prediction of the total system end-to-end 
response J 2 e HB+P (u;) using the double- well protein model. The (a) folded and (b) unfolded states are analyzed separately, with 
the simulation results plotted as symbols and the theory as solid lines. The latter uses the true values of fip = 0.03 fj,o in the 
folded state and /xp = 0.12 no in the unfolded state. Blue/red results denote the Re/Im parts respectively. For comparison, the 
protein response J^(u) (dashed line) is also shown for each well, based on the local Fokker-Planck numerical solution. Inset 
to (b): for the unfolded state, a zoomed-in section of the ImJ™ +p (w) peak, with simulation (symbols) and theoretical (red 
curve) results. The cyan curves are theoretical results with /ip different from the true value: from left to right, ^ip = 0.04, 0.08, 
0.016, 0.24 no- ' 



To implement this approach, we need to extract "single-well" time series from the full simulation data, a portion 
of which is shown in Fig. IS.2T d) . Since what we measure is the end-to-end behavior of the total system, we choose 
the separation line between wells to lie at the barrier position in the full system free energy profile U(z^ B ~ p ) = 
— ksT log T'(z^ B+p ) [Fig. IS.2f e)]. Segments of the time series falling on either side of the separation line will be 
assigned to well 1 (folded) or well 2 (unfolded), with the following exceptions: points that are within a certain time 
window ±5t of a transition (i.e. a crossing of the separation line) are excluded. Setting 5t to 10 6 r, longer than the 
relaxation times within the wells, this exclusion ensures that the resulting single-well time series are not significantly 
perturbed by memory effects from the transitions. Note that times both before and after each transition are excluded 
in order to preserve the time-reversal symmetry of the equilibrium time series. For the time series in Fig. IS.2f d). 
the parts assigned to well 1 and 2 are colored red and blue respectively, with the remainder shaded gray. The same 
method can be generalized to a more complicated free energy landscape, for example in a protein with intermediate 
states, to get a single-well time series corresponding to every state. 

The end-to-end distributions from the single- well series yield corresponding energy profiles Ui(z 2BB+F ) and 



U 2 (z, 



2HB+P* 



, shown in Fig. IS. 2( f). After static deconvolution with the double-HB distribution from Fig. 2(b), we 



get profiles in terms of the protein end-to-end distance: Ui(zf e ) and U2(z^ c ). Near the minima these are identical 
to the original Up, but have anharmonic corrections as they approach infinity at the separation line. The dynamic 
deconvolution analysis for each well is analogous to the simple parabolic case described in the main text, except 
instead of the single Lorentzian we will use the corrected jf e based on the Fokker-Planck numerical solution with 
Ui(ze e ) and [^(z^). This yields a generalized form in each well: Jf c (uj) — J^ Q . a _." , with the multiple Lorentzians 
accounting for the anharmonic corrections. For the parabolic case, C\ = ksT/k, A± = fxpk, and there were no other 
terms. Given a /j,p for the well, the parameters {C a ,A a } are known from the Fokker-Planck solution. (In fact only 
the A a depend on the diffusion constant; note that J^,(0) = Tl a PC a — (3((z p c ) 2 ) is an equilibrium average, so the 
C a must be independent of /Up). 

As in the single-well example, applying the convolution of jf c with the HB response functions, and comparing with 
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the simulation results, shows excellent agreement [Fig. IS.3| . Here the theoretical curves (red lines) for each well are 
based on setting /ip to the values at the minima. For the unfolded state, the inset in Fig. lS.3l shows how this theoretical 
result shifts if /ip is displaced from the true value of 0.12 /io. The sensitivity of the theoretical fitting allows one to 
numerically extract good estimates for /ip in each well: /ip = 0.033 /io (folded) and /ip = 0.111 /io (unfolded), where 
the exact values are 0.03 /io and 0.12 /io respectively. Given the complications associated with excluding interwell 
hopping (and the nonconstant diffusion profile across the barrier), it is remarkable that we can still fit to within « 10% 
of the true values. If one was interested in getting just an estimate of /ip in each well, without necessarily getting a 
perfect fit for the total system response, the anharmonic corrections could be ignored and the single Lorentzian form 
used instead of the exact J^ c in the fitting. The resulting values for /ip have a comparable accuracy. 



VII. ACCOUNTING FOR EXPERIMENTAL LIMITATIONS: NOISE, DRIFT, AND FINITE TIME 

RESOLUTION 

Practical implementation of the deconvolution theory on experimental time series must take into account possible 
errors due to instrumental limitations. We will focus on the optical tweezer apparatus in particular, since this is the 
example analyzed in the main text. However the techniques discussed below are applicable in a generic experimental 
setup. For optical tweezers, effects like drift and noise which enter into the measured output are typically divided into 
"experimental" and "Brownian" categories [13]: the former arises from the instrumental components, while the latter 
is related to thermal fluctuations of the beads and the objects under study. Equilibrium Brownian noise is completely 
captured within the deconvolution theory: the linear response formalism provides a quantitative prediction of exactly 
how each thermally fluctuating part of the system will contribute to the total signal. We will thus concentrate on 
experimental effects, and outline how they enter into the theoretical deconvolution procedure. 

For a perfect instrument, the raw data collected by the experimentalist would be a "true" time series z a (ti) at 
discrete time steps ti, where a refers to either one of the bead positions (a = B,L or B,R), or the bead-bead 
separation (a = ee). The sampling time t s = U + i — U is the time resolution of the apparatus. In any realistic 
scenario, the actual measured time series is z a (ti) = z a (ti) + r] a (ti), where r] a (ti) is an extraneous signal due to the 
instrument. We will consider two possible contributions to rj a (t): a) white noise, or any high-frequency random signal 
that is uncorrelated at time scales above t s ; b) a low- frequency signal related to instrumental drift over an extended 
measurement period. The latter can result from environmental factors like slow temperature changes and air currents 
that affect the laser beam position. As a side note, one of the advantages of the dual-trap setup illustrated in the 
text is a dramatic reduction of such drift effects, relative to earlier designs involving a single trap and biomolecules 
tethered to a surface [13J. Whereas the beam position of a single trap will move over time relative to the surface, the 
dual traps are created from a single laser, and hence any overall drift in the beam position will not affect the trap 
separation. (Though smaller issues may exist, like the effective mobility of the beads shifting as the distance from 
the sample chamber surface varies.) In any case, the low drift of the dual-trap setup allows for long, reliable data 
collection periods (i.e. on the order of a minute in the leucine zipper study of Ref. |14|). 

The final experimental effect we will describe in this section relates to the time resolution: the sampling period t s 
cannot be made smaller than a certain value, dependent on the instrumentation. Not only is the frequency of data 
collection limited, but because the electronics have a finite response time, there will also effectively be some kind of 
averaging of the true signal across the t s time window between measurements. Both of these effects can be incorporated 
into the deconvolution analysis, and we will illustrate this concretely through the toy model protein simulation data 
discussed in the main text. These simulations had data collection at intervals of 0.3t =1.2 ns, in order to clearly 
illustrate that the theory is valid over a wide frequency regime encompassing the characteristic fluctuations of all the 
system components. However we can mimic the experimental apparatus by resampling and averaging the data over 
windows of size t s — 0.01 ms (the resolution of the Gebhardt et. al. [14[ setup). The deconvolution procedure works 
for the reduced frequency range (up to lo s = t^ 1 ), and the main quantity of interest (the protein diffusivity) can 
still be extracted with good accuracy. The reason for this is that the protein characteristic frequency in our example 
(after being slowed down by the handles and beads) falls under the cutoff ui s . In general, any component fluctuation 
modes with frequencies oj < w s can be deconvolved using the theory, even though we lack information about higher 
frequency modes above the cutoff. Though we may not be able to directly observe the modes of certain objects — like 
the short DNA handles whose characteristic frequencies are ~ (9(1 ns) — this does not impede us from exploiting the 
full physical content of the time series in the frequency range below uj s . 

We begin our detailed discussion with the white noise and drift effects that contribute to rj a (t). 
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A. White noise and drift 

Consider the situation where our time series is contaminated by a Gaussian white noise signal rj a (ti), with zero 
mean and correlation {rj a (ti)r] a (tj)) = 7^. We will also assume that the white noise is to a good approximation 
not correlated with the true component signal: {rj a (tj)z a (tj)) = 0, Vi,j. From the measured time series z a {ti) — 
z a {U) + T] a {U) the main quantity which enters the deconvolution theory is the MSD function: 

A a (U) = ((z a (U) - ~z a (0)) 2 ) = ((z a (U) - z a (0)) 2 ) + (( Va (U) - Va (0)) 2 ) 

= A a (ti)+27, 

where A a (ij) is the true MSD value. Thus the white noise induces a uniform upward shift of the MSD, which is 
irrelevant to the theoretical analysis, since we are interested only in the MSD slope, which determines the linear 
response function J a (t) — ((3/2)dA a (t)/dt — J a (t). However this assumes that we can collect an infinite time series, 
reducing the statistical uncertainty of our expectation values to zero. In reality, our sampling is carried out only N 
times, over a finite time interval T s — Nt s . Our calculated A a (ti) for this time interval will differ from the T s — > 00 
result by some random amount characterized by a standard deviation a(A a (ti)). The standard error analysis for 
correlation functions [f 5, 16] (which assumes Gaussian-distributed fluctuations) yields an approximate magnitude for 
a(A a (ti)), up to lowest order in 7: 



<A a (U)) = Ja^ 



lets 7 , Q (t_ 



AT a 1 a V J- a 



(S.35) 



Here A ~ O(l) is a numerical prefactor (which may depend on £,-), T a — (z 2 ) — (z Q ) 2 , and r Q = 
2T~ 2 L dt ((z a (t)z a (0)) — (z a ) 2 ) 2 is an effective correlation time. Since A a (t) can be modeled as a sum of terms 
exponentially converging toward the long-time limit A Q (oo) = 2T a , the value of r a is on the order of the longest 
decay time. Eq. (|S.35[) allows us to see how to correct the effects of white noise: if in the absence of noise (7 = 0) 
we could achieve a certain statistical uncertainty by taking Nq steps, the same uncertainty could be obtained in the 
presence of white noise by taking a larger number of steps, N « A*o(l + 16£ s 7/AraT a ). 

Generally the larger the correlation time r Q for the specific MSD of interest, the longer the sampling time necessary 
to achieve a certain level of precision (i.e. a self MSD, involving the motion of the individual trapped beads, will 
typically have a slower relaxation time than an end-to-end MSD). If the largest r a is ~ 0.1 ms (an upper bound 
estimate for the dual-trap setup discussed in the main text), t s = 0.01 ms, and T a ~ 10 nm 2 is the characteristic scale 
of the MSD function, then a 1% precision, <j(A a (ti))/T a = 0.01, requires a sampling time of T s ~ 1 s (with no white 
noise). The same precision with a noise strength of 7 = 1 nm 2 would need ~ 16% more sampling time. 

In contrast to the high frequency noise modeled above, imagine that we have a low frequency contamination induced 
by slow drift: rj{ti) = ti/td + t/q for some constants t^ and 770. The measured MSD becomes: 



A a (ti) = ((z a (ti) ~ S a {0)) 2 ) = ((z a (U) - z a (0)) 2 ) + ((va(U) ~ V a (0)) 2 ) 

= A a (ti) + (U/t d ) 2 . 



, (S.36) 



Thus instead of reaching a plateau at large times, the MSD continues to increase oc t 2 . If any deviation of this kind 
is observed in the experimental time series, the simplest solution is to fit the long-time MSD, extract the functional 
form (i.e. the slope l/t^), and subtract the drift contribution from A a (ti) to recover A Q (tj). Alternatively, since the 
drift time scale is typically larger than any characteristic relaxation time in the system, t^ S> t q , one could collect 
data from many short runs of length T s , where r a <C T s -C t c \ . 

B. Limited time resolution 

In addition to the possibilities of noise and drift modifying the signal, we have to consider that the measuring 
apparatus will have some finite time resolution t s , defined as the interval between consecutive recordings of the data. 
The measured value, z a (ti), even in the absence of noise contributions r}(ti), only approximately corresponds to the 
instantaneous true value Z a (ti). Because of the finite response time of the equipment, it is more realistic to model 
z a (ti) as some weighted average from the previous t s interval: 

z a {ti) = W- 1 dsz a {U + s)w{s), W= dsw(s), (S.37) 
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for some weighting function w(s). In the present discussion we will ignore any contributions from white noise and drift 
in the time-averaged signal, since these can be corrected for in the same manner as described above. The measured 
autocorrelation function R a (U — tj) = (z a (ti)z a (tj)) between two time points ti and tj is related to the true function 
Ra(U — tj) = (Za(ti)z a (tj)) as follows: 



/o M 

ds / ds' R a (t t + $-tj- s')w(s)w(s') 
-t B J -t. 



/oo 
dt' B(t')R a (U -tj -*'), 
-OO 



(S.38) 



where we introduce the combined weighting function: 

1 „,. . . s~, ... x r lM ,., ft' ~t\ ft 1 + t 

|*|—2t 



B(t) = 2VJF2 e (* + '<>) e H + *■) / _ dt ' w ( — ) w ( ~2~ ) ' ( S ' 39 ) 



While the precise weighting function may vary depending on the apparatus, a reasonable approximation is that 
w(s) sa 1, or that the signal is directly an average over the t s interval. Plugging this into Eq. (|S.39[) . with W — t s , we 
find: 

B(t) = e(t + t s )@(-t + u)j ( 1 -t\ (s - 40) 

The autocorrelation R a (t) enters the deconvolution analysis through the MSD A a (i) = 2(R a (0) — R a (t)), and 
its derivative J a (t) = (/3/2)dA a (t)/dt = —j3dR a {t)/dt. If R a (t) is expressed as a sum of decaying exponentials, 
Ra(t) = J2i Ai exp(~ Ait) , then this corresponds to J a (t) = ^ C* ex P( — ^ii) with the coefficients d = f3AiAi. 
Assuming this underlying exponential form, the relationship between the measured and actual R a (t) can be calculated 
using Eqs. (fST38T) and fl£30j}: 

R a ( t ) = r df B(f)R a (t - f ) = e 2 ( cosh (^ s )-i) Aie _ Ait s ^ Aie _ Ait ^ t ^ ts (s 41) 

J-oo A ih . 

Thus the only modification is in the coefficient of each exponential term, At = p(Ait s )Ai. which gets a prefactor 
p(x) = 2(cosh(a;) — l)/x 2 > 1 dependent oni = Ajt s . The restriction to times t > t s comes from the fact that R a (t) 
cannot be measured for times smaller than the experimental resolution. The prefactor p(x) is only appreciably larger 
than 1 for relaxation times A" 1 comparable to or smaller than t s . However, even for A,^ 1 = t s , which is roughly the 
smallest relaxation time we can realistically fit from the measured data, p(l) — 1.09, so the averaging effect is modest. 

Thus when we fit a sum of exponentials to the J a (t) functions derived from the experimental time series, we are 
effectively calculating Ci = (3AiAi and Aj. To recover the actual d, we just divide out the prefactor: Cj = Ci/p(Ait s ). 
In this way we correct for the distortion due to time averaging. 

To illustrate the effects of the experimental time resolution, and how the above corrections would work in practice, 
we will redo the deconvolution analysis for the toy protein example described in the main text. However, instead 
of using the original simulation time series, whose data collection interval was 1.2 ns, we will average the data over 
windows of width t s — 2500r = 0.01 ms, and use the time series of these mean values, spaced at intervals of t s (with 
a total sampling time T s fa 0.5 s). We thus roughly mimic the apparatus of Ref. [14]. Fig. IS.4f a) shows time-domain 
response functions derived from the coarse-grained data: J s 2 f f B (£), from the first deconvolution step, involving the 
beads and handles alone; ^oc HB+P W> from the second step, where the toy protein is included in the system. The 
fits, drawn as solid lines, involve only single exponentials, since these are sufficient to capture the behaviors in the 
restricted time range t > t s . From the Fourier transform of the fit results for J 2 ™ P (t) (with the Ci — > Ci correction) 
we can get the true frequency space linear response of the 2HB+P system, J™ B+p (o;). The real and imaginary parts 
of this response are drawn as thick solid lines in Fig. IS.4f b) and (c) respectively. As expected, we can calculate this 
response only up to the cutoff frequency u> s = t^ 1 = 4 x 10 _4 r _1 of the measuring equipment. However, it coincides 
very well with the theoretical prediction (thin solid lines, taken from Fig. 3 of the main text). 

Even though we are limited to the range u < u) s , the convolution rules still work: a relation like Eq. [3] in the main 
text is valid at all individual values of ui. Though we see only a restricted portion of the total system end-to-end 
response J™ (o>), it does exhibit non-trivial structure for u < w s , i.e. the downturn in the real part and leveling 
off of the imaginary part at u> > 10 _4 r _1 . Deconvolution allows us to extract the physical content of this structure: 
when numerical fitting is carried out over the range ui < w s , we find toy protein parameters: /ip = 0.0472 ± 0.005 /io, 
fc P = 0.0201 ± 0.0001 k B T/a 2 , within 6% of the exact values (0.05 // and 0.02 k B T/a 2 ). 
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FIG. S.4. (a) Time-domain response functions J s 2 e ™(£) (circles) and J™ B+P (t) (squares) calculated from coarse-grained sim- 
ulation time series of the 2HB and 2HB+P systems respectively. The coarse-graining involved averaging over intervals of 
i s = 2500r = 0.01 ms, and using the time series of mean values. The solid lines are the results of single-exponential fitting to 
the data. The original fit results based on the fine-grained time series, J^p(t) and J™ B+F '(£), are shown as dashed lines for 
comparison. (b,c) The thin dashed and solid lines are the theoretical predictions for the Fourier-domain end-to-end responses 
Je e {uj) and J™ b+p (cj) respectively; panel (b) shows the real part, panel (c) the imaginary part. These are taken from Fig. 3 in 
the main text. The thick solid lines are Re and Im J™ B+P (aj) as determined from the time-domain best-fit in panel (a), after 
correcting for the averaging effect. The vertical dotted line marks the frequency cutoff ui B — t^ 1 for the coarse-grained data. 

What is the main reason behind the relative success of the fitting in extracting the protein parameters, despite 
the limited time resolution? Fig. IS.4f b.c) also shows the theoretical response J^ e (uj) for the protein alone, which 
has a characteristic frequency /ipfcp = 10 _3 t _1 , above the equipment cutoff w s = 4 x 10 r _1 . (Graphically, the 
peak in J^JjS) falls to the right of the vertical dotted line marking ui s .) However, the protein dynamics is modified 
by the presence of the handles and beads in the experimental system, in a manner quantitatively described by our 
convolution theory. As a result, the characteristic frequency of the protein within the experimental apparatus (i.e. 
the peak position in the imaginary part) is shifted to w = 3.6 x 10~ 4 r -1 , just within the cutoff. Hence we are able 
to fit for the protein properties. The relative accuracy of the fitting is notable, since we are essentially at the very 
limit of resolvability. This example illustrates a general principle: whatever dynamical features a protein (or any 
other component) exhibits at uj < lo s within the full experimental setup, our deconvolution theory should be able 
to extract. For this purpose, knowledge of the system behavior for uj > uj s is not necessary. In the current case, 
the theoretical 2HB+P end-to-end response has a contribution at u > 10 _2 r _1 due to the DNA handle fluctuation 
modes, i.e. the shoulder seen in Fig. IS.4f b.c). These modes are experimentally inaccessible at our time resolution t s , 
so the coarse-grained time series reveals nothing about their properties. However, this does not stop us from applying 
the theory at uj < lo s . 
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